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Abstract 

In this paper, we consider a quantum model of gravitation interacting with a Born-Infeld(B-I) type 

scalar field (p. The corresponding Wheeler-Dewitt equation can be solved analytically for both very 

large and small cosmological scale factor. In the condition that small cosmological scale factor tend 

to limit, the wave function of the universe can be obtained by applying the methods developed by 

Vilenkin, Hartle and Hawking. Both Vilenkin's and Hartlc-Hawking's wave function predicts nonzero 
Oh! 

cosmological constant. The Vilenkin's wave function predicts a universe with a cosmological constant 
as large as possible, while the Hartle-Hawking's wave function predicts a universe with positive cos- 
mological constant, which equals to j. It is different from Coleman's result that cosmological constant 
is zero, and also different from Hawking's prediction of zero cosmological constant in quantum cos- 
mology with linear scalar field. We suggest that dark energy in the universe might result from the 
B-I type scalar field with potential and the universe can undergo a phase of accelerating expansion. 
The equation of state parameter lies in the range of —1 <w< — |. When the potential V((p) = j, our 
Lagrangian describes the Chaplygin gas. In order to give a explanation to the observational results of 
state parameter w< — 1, we also investigate the phantom model that posses negative kinetic energy. 
We find that weak and strong energy conditions are violated for phantom B-I type scalar field. At 
last, we study a specific potential with the form Vb(l + vo in phantom B-I scalar field in detail. 

The attractor property of the system is shown by numerical analysis. 

Keywords: Born-Infeld type scalar field; Quantum cosmology; Dark energy; Phantom field 
PACS: 98.80.Cq, 04.50. +h 
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I.INTRODUCTION 



Today, one of the central questions in cosmology is the origin of the dark energy. Many 
candidates for dark energy have been proposed. Among these models, the important ones are 
cosmological constant, " quintessence" [7] , "K-essence" [8] and "tachyon field" [9]. On the other 
hand, nonlinear Born-Infeld theory has been considered in string theory and cosmology 

In this paper, we consider the quantum cosmology and dark energy model of nonlinear 
Born-Infeld type scalar field. The corresponding Lagrangian has been first proposed by 
Heisenbergfl] to describe the process of meson multiple production connected with strong 
field regime, as a generalization of the B-I electromagnetic field Lagrangian L BI = b 2 [l — 
\A + i^) F ^ F ^]i 2 ]- The Lagrangian of B-I type scalar field is 

L s = j{l-^Jl- *<P,v<P,u9»u) ~ V(<p) (1) 

Eq.(l) posses some interesting characteristics [3] that nonsingular scalar field solution can be 
generated, and shock waves don't develop under smooth and continuous initial conditions. 
Second when g^ip^cj)^ <C j, by Taylor expansion, Eq.(l) approximates to 

L s = \f v vAv ~ V &) (2) 

the linear theory is recovered. Third, when potential V = j, our Lagrangian represents a 
universe filled with chaplygin gas which maybe unify dark energy and dark matter [5]. 

We consider quantum creation of universe based on the Wheeler-Dewitt(WD) equation 
Hip = in the superspace. This quantum approach applying to cosmology may help us 
avoid the cosmology singularity problem and understand what determined the initial state 
of the universe. In section II, we consider quantum cosmology with B-I type scaler field 
and constant potential V(<p) which describes cosmological constant effectively. We also find 
the wave function of universe by applying the methods developed by Vilenkin and Hartle- 
Hawking[4]. Both Vilenkin's and Hartle-Hawking's wave function predicts nonzero positive 
cosmological constant. It is different from zero cosmological constant that has been predicted 
by Coleman[10] and Hawking [4]. In fact the H-H's wave function with B-I type scalar field 
predicts a universe filled with Chaplygin gas. In section III, we investigate the cosmology 
with B-I type scalar field as dark energy. Section IV is conclusion. 

II. QUANTUM COSMOLOGY WITH B-I TYPE SCALAR FIELD 
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In order to find the solution of the WD equation, we shall apply the minisuperspace 
model- a closed Robertson- Walker (R-W) metric. In the minisuperspace there are only two 
degrees of freedom: the scale factor a(t) and scalar field ip(t). Using Eq.(l) and by integrating 
with respect to space-components, the action S = J i&fgyf—gd^x + J L s ^/^gd 4 x becomes 

S = j - a 2 )adt + j 2vr 2 a 3 [i(l - yjl - Xtp 2 ) - V(<p)]dt = J C g dt + j C s dt (3) 

where the upper-dot means the derivative with respect the time t. From the Euler-Lagrange 
equation 



we can obtain 



dt y d<p> dp {) 



= 7^w (5) 



where c is integral constant. From the above equation we know that cosmological scale factor 
a is very large or small when ip is very small or large respectively. The critical kinetic energy 

\{ < P)'max ls J\- 

To quantize the model, we first find out the canonical momenta P a = dC g /dd = — (S/2G)aa, 



= dC^/dip = 2ir 2 a 3 (p/\/l — Xip 2 and the Hamiltonian H = P a d + P^ip — C g — C s . H can 
be written as the follows 



G „, 3tt r 8ttG 2Tr/ .. 2vr 2 a 3 r XP 2 

H = -^ p « - 4G° [1 - — ° V( "»i -— I'-r^i (6) 

For ip 2 <C t, the Hamiltonian Eq.(6) can be simplified by using the Taylor expansion, and 
the terms smaller than 6 can be ignored, so the Hamiltonian becomes 

tt G n 2 37T 8irG 2 . P% XP 4 

H = -3^ P « " AG a[l ~ -3-^^)] + 4^3-64^ (7) 

If p{X(p 2 ~ 1) is very large, Eq.(6) becomes 

G n 37r . 8ttG n, T ,, . 1,, 
H = -^Ta P °-4G^-— ^)-J» (8) 

The WD equation is obtained from Hip = 0, Eqs.(7) and (8) by replacing P a — > —i{d/da) 
and — > —i(d/dip). Then we obtain 

a 2 p a 1 a 2 A d 4 

W + a da ~ a 2 d^ 2 ~ 16vr 4 a 8 ^ 

and 



U(a,®)]ip = (9) 



» <) . ... . . 



where <1> 2 = AnGif 2 /3 and the parameter p represent the ambiguity in the ordering of factor 
a and d/da in the first term of Eqs.(7) and (8). We have also denoted 

%*) = (|)V[l-^aM$)] (11) 

u{a^) = {^?a\l-^[Vm-\]} (12) 

Eqs.(9) and (10) are the WD equations corresponding to the action (3) in the case of small 
and large ifi respectively. 

Now we take the ambiguity of the ordering factor p = — 1 and set the transformation 
(o/ao) 2 = a, with ao being the Plank's length. Taking the Plank constant h = 1 and the 



speed of light c = 1, ao ~ ^4G/3tt, we obtain from Eq.(9) 

d 2 ip 1 d 2 tp A d 4 V> 

'da 2 ~ a^d¥ ~ 16^a^ d¥ ~ 



U^ = (13) 



where U = (3n / 4G) 2 o,q(1 — aa^a) and a = (8irG/3)Y. Denoting m = ( 16 ^4 a e ) and assuming 
ip(a, <&) ~ Q(a)e~ K ® with K being an arbitrary constant, we can take the Eq.(13) as follows 

d 2 Q .K 2 mK 4 n 

If a(or a) is large, aa 2 = {8ttGV ((p) /3)a 2 ) a 3> 1 and Eq.(14) approximates to 

^+fU%aQ = (15) 

where // = (3ir/4G)a. The general solution of Eq.(15) can be expressed in terms of Bessel 
function and is given by 

Q(.) = ^|(^fV /2 ) (16) 

Q{<j) is an oscillatory function. 

In the next step, we consider the solution of the WD Eq.(10) with ip being large. We still 
keep the ambiguity of the ordering of factor p = —1. By the same transformation (a/ao) 2 = cr, 
Eq.10) becomes 

^-(%K(l-H 2 a 2 *W = (17) 

where H 2 = (8irG/3)[V(<p) - j]. When H 2 a 2 <C 1, Eq.(17) approximates to 

d 2 ib , 3ir , 9 a , , . 

Solving Eq.(18), we get 

iP = Ne-^ a « u = Ne~^ a2 (19) 



where N is a constant. 

We can see that solution(19) is consist with Vilenkin's tunnelling wave function Eq.(29). 
When a is larger, H 2 a 2 = U 2 a^o » 1, and Eq.(17) approximates to 

^ + {%?<%H^ = (20) 

Its solution is 

Y> = V^(^U) (2i) 

where H = (fg)#. The wave function ip given by Eq.(21) is an oscillatory function. 

Next we will consider the cosmology in case of very large (^(correspondingly very small 
a(t)) by Vilenkin's quantum tunnelling approach. Eq.(10) has the form of a one-dimensional 
Schrddinger equation for a "particle" described by a coordinate a(t), which is zero energy 
and moves in a potential u. The classically allowed region is u < or a > H -1 , with H 
being defined in Eq.(17). In this region, disregarding the pre-exponential factor, the WKB 
solutions of Eq.(10) are 

V4 1} (a) = exp{[±i f P{a')da'} T ^} (22) 
The under-barrier(a < H -1 , classically forbidden or Euclidean region) solutions are 

V4 2) (a) = exp{[± [ H \P(a')\da'}} (23) 

J a 



where P(a) = y/\ — u(a)\. 

The classical momentum conjugate to a is P a = —ad. For a > H -1 , we have 

(-^)#(a) = ±P(a)#(a) (24) 

and thus ip^\a) and ip+\a) describe the expanding and contracting universe respectively. 
The tunneling boundary condition requires that only the expanding component should be 
present at large a, 

iP T (a > H- 1 ) = V^V) (25) 
The under-barrier wave function is found from WKB connection formula 

Ma<H- 1 )=^\a)- l -^\a) (26) 

(2) (2) 

The increasing exponential ip_ (a) and the decreasing exponential ip + (a) have compara- 
ble amplitudes at the nucleation point a = H^ 1 , but away from that point the decreasing 



exponential dominates 



Ma < H' 1 ) « ^){a) = exp[^(l - H 2 a 2 )i] (27) 

The " tunneling amplitude" (probability distribution for the initial values of V in nucleating 
universe) is 

V't(O) 

From Eq.(28) we obtain the result that the tunneling wave function predicts a nucleating 
universe with the vacuum energy(i.e cosmological constant) as large as possible and critical 
kinetic energy as small as possible (defining the kinetic energy when a = as critical kinetic 
energy, which equals from Eq.(5)). If H 2 a 2 <C 1, by Taylor expasion Eq.(27) becomes 

, . „ 7T 37T0 2 . 

Ma) ~ e ^[^2 " "4G"] ( 29 ) 

Comparing Eqs.(19) and (29), we find that the only difference is just an unimportant pre- 
exponential factor. 

The Hartle-Hawking(H-H) no boundary wave function is given by the path integral 

$hh = J[dg][d^ SE ^ (30) 

In order to determine ipHH, we assume that the dominant contribution to the path integral 
is given by the stationary points of the action(the wormhole instantons) and evaluates ipHH 
simply as ipHH oc e~ SE ^ saddie -p° int . When (0) 2 ~ j, from action(3) we can obtain 

S = J H [(1 - a 2 )adt - J 27r 2 a 3 H 2 dt] (31) 

where H 2 = ^-(v — j). The corresponding Euclidean action Se = —i(S) CO ntinue is 

S E = j-^[l + ^) 2 ]adr + J 2n 2 a 3 H 2 dT (32) 

where r = it. From action(31), we can obtain that the a(t) satisfies the following classical 
equation of motion 

_(^)2_ 1 + F 2 a 2 = () (33) 

The solution of Eq.(33) is the de sitter space with a(t) = H~ l cosh{Ht). The corresponding 
Euclidean version (replacing t — ► — it) of Eq.(33) is 

{ ^Lf_ 1 + H 2 a 2 = Q (34) 

ar 



The solution of Eq.(34) is 



o(t) = H^siniHr) 



(35) 



Using Eqs.(32)(35), we obtain 



iPhh(o) oc exp[- 



TT 



(l-ffV)i] 



(36) 
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The only difference between the H-H's wave function(36) and Vilenkin's wave function(27) 
is the sign of the exponential factor. The H-H's wave function(36) gives the probability 
distribution 



The H-H's distribution(37) is the same as Vilenkin's one(28), except a sign of the exponential 
factor. The distribution(37) is peaked at V — (1/A) = and it predicts a universe with a 
positive cosmological constant j. This is a very significative result, it is different from zero 
cosmological constant predicted by Coleman and Hawking. 

When V = j for our lagrangian, it represents a universe filled with Chaplygin gas [5]. In 
fact H-H's wave function predicts a unverse filled with Chaplygin gas. Next, we will discuss 
dark energy model of B-I type scalar field and the relationship between our Lagrangian and 
Chaplygin gas. 

III.Dark Energy Model of B-I Type Scalar Field 
1. The Model With Lagrangian ±[1 - y/1 - Xg^p^p^\ - V 

In the spatially flat R-W metric, Einstein equation G^ v = KT^ U can be written as 



7T 



(37) 




(38) 



2 ^ + (it)2 = KTl = KTl = KT$ = -Kp 



(39) 



Substituting Eq.(38) into Eq.(39), we get 



(40) 



where 



The energy density p s = Tq and pressure p s = — T/ = — T 2 2 = — T| are defined as following 



T . 1 , 1 5 + VI + Ac 2 a- 6 

Ps = y -A + vr=^ = a (42) 

p s = i[l-^V]-7 (43) 

where the upper index "." denotes the derivative with respect to t, and V(tp) is taken Vq as 
a constant. We define B = XVq — 1. When a(i) = the kinetic energy ^ = ^ is critical 
value from Eq.(5). According to Eqs.(38) and (42) we get 

= ^l + \ c 2 a ^ + B } (44) 

Ob OA 

When a(t) is very small(a(i) — ► 0), Eq.(44) approximates to 



& = J£- (45) 
V 3A2a 



2 



a oc t3 (46) 
When a(i) increases little by little until a 3> (Ac 2 )e, Einstein equation(44) becomes 



From Eq.(47) we obtain 



(^ = f[Vo + ^-] (47) 



.. KV Kc 2 _ 5 

a = -^-a - —a 5 (48) 

„2 , 1 



From the above equation we find that when a > (f?-) 6 , the universe is undergoing a acceler- 
ated phase. Integrating Eq.(47) we obtain 



/ r 2 ( p y/12KV t _ 11 

qSpc /_£_^ e (49) 



When a — > oo, Eq.(44) becomes 



(50) 



From the above equation, we see that the universe is undergoing a inflation phase. However, 
the dark energy model interpolates between a dust dominated phase in the past and a de-sitter 
phase at the late time. 

Especially when Vq = j, our Lagrangian describes the Chaplygin gas that was proposed 
as a model for both dark energy and dark matter in the present universe[5]. Correspondingly, 
the density and pressure are 

VI + Ac 2 a- 6 



Substituting Eq.(51) into Eq.(38), we obtain the solution of Einstein equation(38) 

t = ^[ln Y + - K + ^ -2 a rctan(Vl + c2Aa^)i] (53) 
6 (^ + T«" 6 ) 3 -(i)" 



When a(t) is very small(a(i) — > 0), the density is approximated by 

c 



a 3 A 2 



(54) 



that corresponds to a universe dominated by dust-like matter. For large values of the cosmo- 
logical scalar factor a(t) it follows that 

1 . . 

Pc ~ j (55) 

p c ~-i (56) 

which in turn corresponds to a universe with a cosmological constant j(i.e., a de-sitter uni- 
verse). At early time, i.e., the cosmological scale factor a(t) is small, p c ~ — S-r, which 

a 3 A2 

corresponds to a dust like dominated universe. At late time, i.e., the cosmological scale fac- 
tor is large, p c ~ t, which corresponds to a cosmological constant like dominated universe. 
Therefore the Chaplygin gas interpolates between a dust dominated phase in the past and a 
de-Sitter phase at late time. 

When potential is taken to zero, the density and pressure are 



A [l-Vl"V] (57) 

Ps = \[l-\]l-^ 2 ] (58) 
Substituting Eq.(5) to the above two expressions, we have 

U7= — = = (59) 

Ps Va 6 + Ac 2 V ' 

and can see 

< w < 1 (60) 

So, there is no accelerated expansion in the universe of B-I type scalar field without potential. 
For B-I type scalar field with potential, we have 
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When potential is greater than j and the kinetic energy of (p field evolves to region of (p 2 < ^ , 
p + 3p < 0. The universe undergoes a phase of accelerating expansion. 
We also get 

p+p= tSw >0 (62) 



and 



w = - > -1 (63) 



However, some analysis to the observation date hold that the range of state parameter lies 
in —1.32 <w< — 0.82[6]. From Eq.(63) we know that parameter of state equation is larger 
than -1. In order to give a favor explanation to the observation results, we investigate the 
phantom field model that possess negative kinetic energy and can realize w< —1 in their 
evolution. Next we discuss phantom field model. 

2. The Model With Lagrangian ±[1 - y/1 + Xg^ip,^] - V(ip) 

We consider the case that the kinetic energy terms is negative. The energy-momentum 
tensor is 

T}} = C^Mk - 5£L (64) 



From Eq.(64), we have 



p= _ T ;_-T| = -T| = i-^±^:-F (66) 



Based on Eq.(65) and (66), we can obtain 



p+ " = -tSw (67) 



It is clear that the static equation w< —1 is completely decided by Eq.(67). We can also get 

' +!f TI^-iW" i1 ' (68) 

It is obvious that p + 3p < 0. Eq.(68) shows that the universe is undergoing a phase of 
accelerated expansion. The model of phantom B-I type scalar field without potential V(tp) 
is hard to understand. In this model we can always find p = — , 1 — 4- < and (-) 2 < 0. 
It is unreasonable apparently. In the model of phantom B-I type scalar field with potential 
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V(<p), if V(ip) > j — ^-j===, p is greater than zero. 



First, we consider the case of a specific simple example V = uq = const and no — j 



A -{A > 0). So, Eq.(63) becomes 



Xy/1 + \<p 2 A 

Substituting Eq.(69) into Eq.(38), we have 



+■ 4 ( 69 ) 



(~) 2 = ^ L ' . 2 +y] (70) 



We can obtain from the Euler-Lagrange equation (4) 



= 7W=W (71) 



where c is integrate constant. Substituting Eq.(71) into Eq.(70), we get 



a 



^[Vl-Xc^a^ + A] (72) 



From the above equation, we know the universe is nonsingular because the minimum a m in of 
scale factor is (Ac 2 )^. When the universe scale factor approximates to a m j n , Eq.(72) becomes 



h = v fr (73) 



a = e V— * (74) 



When a — > oo, Eq.(72) becomes 



. 1^4 + 1) 

a = V^A^ a (75) 



a = eV 3A ( 76 ) 

In phantom model with constant potential, the universe is always undergoing a phase of 
inflation at very small and lager scale factor. 

Next we study the cosmological evolution by numerical analysis in the phantom model 
with a specific potentialfll] 

V(<p) = Vb(l + — )e ( "^ } (77) 
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As we consider the phantom field becomes dominant , we can neglect the nonrelativitic 
and relativistic components (matter and radiation) in the universe, then from Eqs.(5,38), we 
have 

(p + v{l + A<^3iq^=L== - I + V(<p)] - v'(cp)(l + A0 2 )§ = (78) 

where the overdot denotes the differentiation with respect to t and the prime denotes the 
differentiation with respect to ip. 

To study an numerical computation, it is convenient to introduce two independent vari- 
ables 

X = ip 
Y = ip 



(79) 



then Eq.(78) can be written 



dX \s 

dt ~ 1 



% = V\X){1 + XY 2 )^ - Y(l + XY^) y /3K[ T7 ±- s -j : + V(X)] 
We can obtain this system's critical point from 



(80) 



— — fi 

dt — u 

dt u 



(81) 



then its critical point is (X c ,0), where the critical value X c is determined by V'{X C ) = 0. 
Linearizing Eq.(80) around the critical point, we have 



dX =Y 

dt (82) 

<f = V"(X C )(X - X c ) - ^(3KV(X C )Y 



the types of the critical point are determined by the eigenequation of system 

e 2 + as + (3 = (83) 

where a = ^3KV(X C ),P = -V"(X C ), the two eigenvalues are e x = -V^V(x c )+^3KV(x c )+w"(x c 
e 2 = -V3KV(x c )-^3KV(x c )+4V"(x ])__ ^ & pogitive potentialS; jf V"(X C ) < 0, then the 
critical point (X c ,0) is a stable node, which implies that the dynamical system admits at- 
tractor solutions. We can also conclude that if a potential possesses the general properties: 
V(X C ) > 0, V'(X C ) = and V"(X C ) < 0, then our phantom model with this potential will 
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have a attractor solution and predict a late time de-sitter like behavior (w^ = — 1). Substi- 
tuting Eq.(77) into Eq.(80), we obtain 



d2L - v 

dt ~ 1 

Vo 



(84) 



From Eqs(81,84) we obtain the critical X c = and V"(X C ) = < from Eq(77). 
Therefore this model has an attractor solution which corresponds to its attractor regime, the 
equation of state w < — 1. To solve this equations system via the numerical approach, we 
re-scale the quantities as x = s = (i^gA)"^ y = \/XY. Then Eq.(84) becomes 

dx _ „, 
ds ~ 

f s = _ 7se (-x)(l + y2)| _ 0oy( i + y 2 )y3K[^i=-l+7(l+x)e(-)] (85) 

where we set K = 1 and 7 = VoA is parameter. The numerical results with different initial 
condition are plotted in Figs.l — 3 and the parameters (po = \/0T, 7 = 3. 



1.5- 



X 




2 4 6 8 10 



Figl.This plot shows the evolution of the scalar field in difference initial condition, 
solid line is for (fi n = 0.7(po, dotted line is for (fi n = ipo, dashed line and dot-dashed 
line for ip; m = 1.3c^o 5 1-7(A) respectively, they are all plotted for a fixed value of 
Din = 0.1 . 
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2 4 6 8 10 



Fig2. The evolution of w with respect to s, the initial condition is the same as figl. 













0.5- 


/ 










0- 






) 






-0.5- 


\ \ \ \ 
\ \ \ ^^^^w 


J j 




/ 

/ 


/ 

/ 




y -1- 
















./' 






-1.5- 


\ \ 










-2- 












-2.5- 

















0.5 1 

X 




1.5 



Fig3. The attractor property of the system in the phase plane, the initial condition 
is the same as figl. 

As we know, when A — > 0, our model turns to be phantom quintessence model. In order 
to see the nonlinear effect, we plot the phantom quintessence model with our model in fig4 
and fig5. 
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X 




1 2 3 4 5 6 



Fig4. the evolution of scalar field with respect to s, solid line is phantom quintessence 
model, dotted line and dashed line are nonlinear B-I phantom field, dotted line is for 
A = 2/3,dashed line is for A = 1/3. 




-1 .4 - 
-1.5- 
-1.6- 

1 2 3 4 5 6 

s 

Fig5.the evolution of w with respect to s, solid line is phantom quintessence model, 

dotted line and dashed line are nonlinear B-I phantom field, dotted line is for A = 

2/3, dashed line is for A = 1/3. 
From figl,fig2 and fig3, we can easily find the system admits a attractor solution. The 

equation of state parameter w starts from regime of near -1, then quickly evolves to the 

regime of smaller than -1, last, turns back to execute the damped oscillation, and reaches 

to -1 eventually for ever. Due to the unusual physical properties in phantom model, the 

phantom field, releases at the distance from the origin with a small kinetic energy, and moves 
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towards the top of the potential and crosses over to the other side, then turns back to execute 
the damped oscillation around the critical point. After a certain time the motion ceases and 
the phantom field settles on the top of the potential permanently to mimic the de Sitter-like 
behavior (w^ = — 1). Fig4 and fig5 indeed shows that the nonlinear scalar phantom field will 
reduce to phantom quintessence model when A decreases to zero. The nonlinear effect does 
not affect the global attractor behavior but change the evolution of w and scalar field 4> va - 
details. 

IV.CONCLUSION 

It shows that the tachyon can be described by a B-I type Lagrangian resulting from string 
theory. It is clear that our Lagrangian j[l — y/1 — Xg^ u ip^(p^] — V(ip) is equivalent formally 
to the tachyon type Lagrangian — jyjl — g^ u ^,^,u, with a potential j — V(ip) where we 
re-scale the scalar field as <3? = \^<p. The WD equation are solved analytically for both very 
large and small a(t). For the very small cosmological scale factor, the Vilenkin's tunnelling 
wave function Eq.(27) predicts a nucleating universe with the biggest possible vacuum energy 
V (cosmological constant) and the smallest possible critical kinetic energy The only 
difference between the H-H's wave function Eq.(36) and Vilenkin's wave function Eq.(27) is 
the sign of the exponential factor. Then the H-H's probability distribution in Eq.(37) reaches 
to peak at V — j = 0. It predicts a universe with a postive cosmological constant j. It is 
different from zero cosmological constant that has been predicted by Coleman and Hawking. 

The parameter w of state equation lies in the range of <w< 1 for B-I type scalar field 
without potential. The universe of B-I type scalar field with potential can undergo a phase 
of accelerated expansion. The corresponding parameter of state equation lies in the range of 
— 1 <w< — |. This model admits a late time attractor solution that leads to an equation of 
state w= —1. 

In order to give a favor explanation to the observational results that the range of state 
parameter lies in —1.32 <w< —0.82, we investigate the phantom field model that possess 
negative kinetic energy. The weak and strong energy condition are violated for phantom B-I 
type scalar field. The parameter w of state equation lies in the range of w< —1. When the 
potential V(tp) = uq, the universe is nonsingular and stays in de-sitter phase always. When 
we choose a potential as Eq.(77), the evolution behavior of the state parameter w is shown 
in Fig.2. By numerical analysis, we learn that there is an attractor solution in Fig.3. 
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